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We present a preliminary search for D°-D° mixing using the decays D° -» K+ tt tt°, additionally presenting 
Dalitz-plot distributions and a measurement of the branching ratio for this mode. A new tagging technique 
is used to produce the doubly Cabibbo-suppressed Dalitz plot, which in turn is used to motivate the method 
used for the D-mixing search. We analyze 230.4 fb -1 of data collected from the BABAR detector at the PEP-II 
collider. Assuming CP conservation, we find Rm < 0.054% with 95% confidence, and we estimate that the data 
are consistent with no mixing at a 4.5% confidence level. We present D-mixing results both with and without 
the assumption of CP conservation. 



1. Introduction 

Although A and B mixing are well established, D 
mixing has yet to be observed. As this particular 
mixing phenomenon is sensitive to new physics in a 
complementary manner to the A and B systems, it 
is an essential test of the completeness of the Stan- 
dard Model. We present preliminary results of a new 
mixing analysis using the decays D° — > A + 7r~7r° , tak- 
ing into account the resonant structure of the doubly 
Cabibbo-suppressed contributions. This is the first 
search for mixing using either this decay mode or this 
technique. 

G. Burdman and I. Shipsey have written a thorough 
review of D-mixing predictions [1]; accurate predic- 
tions arc difficult to obtain because of significant con- 
tributions from long-distance effects. 

Using 230.4 fb -1 of data collected from the BABAR 
detector at the PEP-II collider, corresponding to ap- 
proximately 300 million cc events, we obtain a pure 
sample of D° candidates by reconstructing the decays 



D* 



D°- 



K^TT (+ CC). 



(1) 



The flavor of the D° candidate at production is carried 
by the charge of the associated irf . The large sam- 
ple of Cabibbo-favored (CF) decays D° -> A^tt+tt 
is used both as a means for reducing systematic 
uncertainties in the probability density functions 
(PDFs) used to describe the wrong-sign (WS) decays 
D —> A + 7T~7r° and as a normalization mode when 
determining the WS decay rate. 

We separate correctly reconstructed decays from 
background, and distinguish doubly Cabibbo- 
suppressed (DCS) contributions from CF mixed 
contributions, by means of an (unbinned) extended 
maximum likelihood fit. PDFs are fit to the three 
distributions {rn^^o, Am, ^^0}, where rriKwK is 
the invariant mass of the D° candidate, Am is the 
invariant mass difference between the D* + and D° 
candidates, and tjc^o is the candidate decay time. 



2. D° -> K+tt-tt Branching Ratio 

In a nonleptonic search for D mixing, DCS contri- 
butions obscure signs of mixing in the final state. To 
the extent that the DCS rate is low for a particular 
mode compared to the corresponding CF rate, there 
is greater sensitivity to a potential mixing signal. The 
branching ratio 



R 



T(D° -> A+7T-7T ) 
T(D° -V A-TT+TT ) 



(2) 



is measured using a maximum likelihood fit to the dis- 
tributions {m^^o, Am}, and we find the preliminary 
result 

R = (0.214 ± 0.008 (stat) ± 0.008 (syst))%. (3) 

This result is consistent with that reported by the 
Belle Collaboration last year [2]. Comparing this to 
the corresponding branching ratio for the decay D° — * 
A+7T- [3], 



A(Att) = (0.362 ± 0.029)%, 



(4) 



we find that an analysis of D° — > K + ir~ir° may in 
fact have more sensitivity to mixing than the standard 
analysis of the decays D° — > K + ir~ . We note that 
while the branching ratio in Eq. 3 may contain contri- 
butions from D mixing (known to be small) , the level 
of these contributions cannot be determined without 
an analysis of the WS decay-time distribution. 



3. Event-Level Tagging and DCS 
Resonance Contributions 



7T° may 



Just as an analysis of the decay D° — * K + tt 
offer greater sensitivity to a mixing signal than that of 
D° — > K + ir~~ because of the lower relative DCS con- 
tribution, regions of the decay phase space (i.e., the 
Dalitz plot) may have greater sensitivity to mixing 
to the extent that the DCS contribution in those re- 
gions is relatively low. However, obtaining an accurate 
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Figure 1: The m KlT ^o distribution with fitted PDFs overlaid for the WS data sample used to search for D mixing (left) 
and for the sample used to obtain the Dalitz-plot distributions (right). The sample on the right requires a second, 
event-level, flavor tag. Both distributions require 0.1444 < Am < 0.1464 GeV/c 2 . 




Figure 2: The Am distribution with fitted PDFs overlaid for the WS data sample used to search for D mixing (left) 
and for the sample used to obtain the Dalitz-plot distributions (right). The sample on the right requires a second, 
event-level, flavor tag. Both distributions require 1.85 < m Knn a < 1.88GeV/c 2 . 
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Figure 3: Dalitz plot and one-dimensional projection onto rn 2 K7r0 for the CF decay (left) and the DCS decay (right). A 
statistical background subtraction has been applied (i.e. candidates have been weighted), and distributions have been 
corrected for efficiency variations as a function of the Dalitz plot. The shaded histograms in the lower plots are 
projections of the Dalitz plots above them, with error bars superimposed. Candidates have been selected requiring the 
event-level tag. 



distribution of events in the Dalitz plot is challeng- 
ing because of the large number of background events 
still remaining after requiring event-selection criteria. 
In particular, a large component of the background is 
from real D° decays having a misassociated . Since 
this background is from real D° — >■ K~ir + Tr° decays, it 
populates the Dalitz plot with structure that tends to 
obscure that of the signal. This peaking background 
is marked "Bad D*+" in Figs. 1 and 2. 

To suppress peaking background from 
D° — > A'~7r + 7T° decays, we use an event-level 
tag in addition to the primary irf 1 flavor tag. This 
tag is determined by searching the event hemisphere 
opposite that in which the D* + candidate is re- 
constructed for any of the following: K^, nf, e T , 
or fj, T . If the charge of one of these candidate is 
consistent with the hypothesis of hadronization and 
decay of a c quark in that event hemisphere, then the 
event is tagged. The signal efficiency after applying 
this second, opposite-side, tag is 46.4%, while the 
peaking-background efficiency is 10.9%. The effect of 



requiring this tag is shown in Figs. 1 and 2. 



After requiring the opposite-side tag, correcting 
for efficiency variations as a function of the Dalitz 
plot, and performing a statistical background sub- 
traction [4] based on a maximum likelihood fit to 
{tuk^o, Am}, we obtain the Dalitz plots in Fig. 3. 
The CF Dalitz plot is qualitatively different from the 
WS Dalitz plot, which is assumed to contain pri- 
marily DCS contributions 1 . While CF decays pro- 
ceed primarily through the resonance D° — * K~p + , 
DCS decays proceed primarily through the resonance 
D a — > K* + ir~ . We use this observation to maximize 
sensitivity to a potential D-mixing signal as described 
in the next Section. 



x We assume that possible contributions from D mixing can- 
not be distinguished by eye in the Dalitz plots shown. 
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Figure 4: Contours of constant Alog£ = 1.15, 3.0 in 
terms of the doubly Cabibbo-suppressed rate and the 
integrated mixing rate. The upward slope of the contour 
indicates negative interference. 



Figure 5: Contours of constant Alog£ = 1.15,3.0 in 
terms of the normalized interference term and the 
integrated mixing rate, for the D° and D° samples 
separately. The hatched regions are physically forbidden. 



4. Preliminary Results of a D-mixing 
Search 

The two mass eigenstates 

\D h2 )=p\D°)±q\D ) 



(5) 



generated by mixing dynamics have different masses 
(mi, 2) and widths (IY2), and we parameterize the 
mixing process with the quantities 



m,2 — mi 

! r 2 + r 1 : 



y 



(6) 



If CP is not violated, then \p/q\ = 1. For a multibody 
WS decay, the time-dependent decay rate, relative to 
a corresponding right-sign (RS) rate, is well approxi- 
mated by 



r M (t) 



— Rd 



ay y R D (Tt) 
< a < 1, 



x' 2 + y' 2 



(rt) 2 (7) 



where the tilde indicates quantities that have been in- 
tegrated over the selected phase-space regions. Here, 
Rd is the integrated DCS branching ratio; y' = 
y cos 8 — x sin 8 and x! = x cos 8 + y sin 8, 8 is an un- 
known integrated strong-phase difference; a is a sup- 
pression factor that accounts for strong-phase varia- 
tion over the region; and T is the average width. The 
time-integrated mixing rate Rm — {i a + y' 2 )/2 = 
(x 2 +y 2 )/2 is independent of decay mode and should 
be consistent among mixing measurements. 



The necessity of the suppression factor a can be 
understood as follows. Suppose the strong-phase dif- 
ference between DCS decay and mixing is such that 
in a particular region of phase space there is positive 
interference, while in another region, there is nega- 
tive interference. If one performs a mixing analy- 
sis considering the decay-time distributions of events 
in both phase-space regions simultaneously, then the 
combined decay-time distribution will have a sup- 
pressed interference effect relative to the mixing rate 
Rm, which is nonnegative and invariant across phase 
space. 

In addition to our search for a CP-conserving mix- 
ing signal, we manifestly permit CP violation by fit- 
ting to the D° -> K+tt-jt and D° -> R-tt+tt sam- 
ples separately. We consider CP violation in the in- 
terference between the DCS channel and mixing, pa- 
rameterized by an integrated weak-phase <f>, as well as 
CP violation in mixing, parameterized by \p/q\. We 
assume CP invariance in both the DCS and CF rates. 
The transformations 



ay 



±1 



(ay 1 cos <j> ± fix' sin (f>) (8) 
±2 



(x 2 + y 2 ) 



{x 2 + y 2 ) 



(9) 



are applied to Eq. 7, using (+) for 
T(D° -> K-n+TT )/T{D -> K-tt+tt") and (-) 
for the charge-conjugate ratio. The parameter (3 is 
analogous to a. 

Because decays with mixing will obey the CF-decay 
resonance structure, and we can qualitatively sec the 
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Table I Preliminary signal-candidate yields determined 
by the two-dimensional fit to the {m Kn7T a , Am} 
distributions for the WS and RS samples. Yields are 
shown for (a) the selected phase-space regions used in 
this analysis and (b) the entire allowed phase-space 
region. Uncertainties are those calculated from the fit, 
and no efficiency corrections have been applied. 



D° Cand. 



D° Cand. 



(a) 



(b) 



WS 
RS 



(3.84 ± 0.36) x 10 2 

(2.518 ±0.006) x 10 5 

WS (7.5 ±0.5) x 10 2 

RS (3.648 ± 0.007) x 10 5 



(3.79 ±0.36) x 10 2 
(2.512 ± 0.006) x 10 s 

(8.1 ±0.5) x 10 2 
(3.646 ± 0.007) x 10 5 



differences between the CF and DCS Dalitz plots, we 
choose regions of the Dalitz plot to analyze in order to 
maximize sensitivity to a potential mixing signal. We 
do this by excluding events with two-body invariant 
masses in the ranges 850 < m(Kir) < 950MeV/c 2 
or 850 < m(Kir°) < 950MeV/c 2 . The signal yields 
before and after applying these selection criteria are 
listed in Table I. 

We search for mixing in the WS decay-time dis- 
tribution by constructing a three-dimensional PDF in 
{mjf„o,Am, tKTTir"}- The WS signal PDF in t Knn o is 
a function based on Eq. 7 convolved with three Gaus- 
sians. We fit the WS PDF to the t K7T7r o distribution 
allowing yields and shape parameters to vary 

The results of the decay-time fit, both with and 
without the assumption of CP conservation, are listed 
in Table II. The statistical uncertainty of a partic- 
ular parameter is obtained by finding its extrema for 
A log£ = 0.5; in finding the extrema, the likelihood is 
kept maximal by refitting the remaining parameters. 
Contours of constant Alog£ = 1.15,3.0, enclosing 
two-dimensional coverage probabilities of 68.3% and 
95.0%, respectively, are shown in Figures 4 and 5. The 
likelihood is maximized when evaluating A log C for a 
particular point in two-dimensions. 

We note that A log £ as a function of the quantity 
sign(ay') x Rm is approximately parabolic. The two- 
sided interval -0.054% < sign(ay') x R M < 0.054% 
contains 95% coverage probability; thus, we quote 
Rm < 0.054% as our upper limit on the integrated 
mixing rate under the assumption of CP conservation. 

A feature of Alog£ in one dimension is that it 
changes behavior near Rm = because the inter- 
ference parameters become unconstrained. There- 
fore, we estimate the consistency of the data with 
no mixing using a frequentist method. Generating 
1000 simulated data sets with no mixing, each with 
58,800 events representing signal and background in 
the quantities {mx^^o, Am, tx^}, wc find 4.5% of 
simulated data sets have a fitted value of Rm greater 



than that in the observed data set. We conclude that 
the observed data are consistent with no mixing with 
4.5% confidence. 



Table II Preliminary mixing results assuming CP 
conservation (D and D° samples are not separated) and 
manifestly permitting CP violation (D° and D° samples 
are fit separately). The first listed uncertainty is 
statistical, the second is systematic. Rd is not reported 
for the latter case because irf efficiencies were not 
studied. 



CP conserved 


CP violation allowed 


Rm (0.023 ±° ° {l ± 0.004)% 
Rd (0.164 +° ° f 2 ± 0.012)% 




(0.010 ±°;°»)% 


ay -0.012 1°;°°® ± 0.002 


ay cos (f> 
fix' sin 4> 


-0.012 t° Z 

n nrn + 002 

U.UUO _ 005 




\p/q\ 


9 9 +1.9 
-1.0 
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